Control and femtosecond time-resolved imaging of torsion in a chiral molecule 
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We study how the combination of long and short laser pulses, can be used to induce torsion in 
an axially chiral biphenyl derivative (3,5-difluoro-3',5'-dibromo-4'-cyanobiphenyl). A long, with 
respect to the molecular rotational periods, elliptically polarized laser pulse produces 3D alignment 
of the molecules, and a linearly polarized short pulse initiates torsion about the stereogenic axis. 
The torsional motion is monitored in real-time by measuring the dihedral angle using femtosecond 
time-resolved Coulomb explosion imaging. Within the first 4 picoseconds, torsion occurs with a 
period of 1.25 picoseconds and an amplitude of 3° in excellent agreement with theoretical calcula- 
tions. At larger times the quantum states of the molecules describing the torsional motion dephase 
and an almost isotropic distribution of the dihedral angle is measured. We demonstrate an orig- 
inal application of covariance analysis of two-dimensional ion images to reveal strong correlations 
between specific ejected ionic fragments from Coulomb explosion. This technique strengthens our 
interpretation of the experimental data. 



I. INTRODUCTION 

The control of molecules and chemical reactions with 
lasers is one of the main goals of femtochemistry [1 . One 
particular aspect of such studies is the use of laser pulses 
to control the transition from one enantiomer of a chi- 
ral molecule to its mirror form. This topic has been 
the subject of a large number of theoretical studies and 
these efforts are motivated by the intriguing prospect of 
light-induced deracemization [2H81, i.e., creation of enan- 
tiomeric excess. Axially chiral molecules have attracted 
special interest because the potential barrier separating 
an enantiomeric pair is often low so that lasers at moder- 
ate intensities may be applied to facilitate efficient cross- 
ing of the barrier without fragmenting the molecule. For 
axially chiral molecules the stereogenic element is an axis 
(connecting two atoms) rather than a point (a single 
atom) and the reaction path separating the two enan- 
tiomers corresponds to torsion around the stereogenic 
axis |9J. 

The first theoretical light-based scheme which suc- 
cessfully showed the production of enanotiomeric excess 
(50.0001% : 49.9999%) involved circularly polarized con- 
tinuous wave monochromatic light [2 . However, a much 
higher conversion rate is needed if the deracemization 
is to be confirmed experimentally - not to mention to 
have practical applications. To this end studies showed 
that almost perfect deracemization could be achieved by 
exposing phosphinothioic acid (H2OPSH) to picosecond 
terahertz laser pulses [3]. In these studies, the key was 
to realize that changes in the torsional coordinate was 
associated with a shift in the permanent dipole moment. 



a feature that was exploited by using polarization shaped 
laser pulses to transfer a racemate into the desired enan- 
tiomeric form. Since this discovery, several different ap- 
proaches have been proposed to control torsional motion 
in axially chiral molecules. These approaches include the 
use of coherent control or phase controlled pulses [10] 
to perform, e.g., laser distillation [4 where excitation 
of one of the enantiomers can be enhanced by coupling 
the ground state to a superposition of two excited states 
with opposite parities, or by cyclic population transfer 
that relies on quantum interference in an effective three- 
level-system [5 . Alternatively, the deracemization can 
be achieved by exploiting that the direction of the elec- 
tronic transition moment vector to an excited state is 
different for the two enantiomers. One way to exploit 
this directionality is a pump- dump scheme [6 , where the 
population is gradually transferred from one enantiomer 
to the other. Another method would be asymmetric ex- 
citation^ where the excited state is repulsive and leads 
to dissociation thereby removing this enantiomeric form 
from the ensemble [7 . A common feature for most of 
these schemes is that they either require or benefit from 
pre-aligned and oriented molecular targets [8 . 

Since alignment and orientation is generally required in 
the enantiomer-control-schemes a more straightforward, 
and thereby experimentally feasible approach, would be 
to rely on alignment concepts alone to control the tor- 
sion. Recently, this has been the focus of several stud- 
ies. These studies fall into two regimes: (i) the adiabatic 
regime where the pump laser that induces torsion is long 
compared to the time scale of torsional dynamics pTHl3] 
and (ii) torsion induced by a short non-adiabatic pump 



pulse [m [15] . For those axially chiral molecules where 
the stereogenic axis and the most polarizable axis co- 
incide it was suggested that torsion could be controlled 
adiabatically by a single elliptically polarized nanosecond 
(ns) laser pulse provided that the torsional barrier is suf- 
ficiently low. This pulse potentially serves to both induce 
3D alignment and torsional motion by interaction of the 
minor axis of the polarization ellipse with the polariz- 
ability components of the two moieties forcing them into 
coplanarity [11, 13 . The feasibility of this approach has 
recently been disputed due to the coupling between over- 
all rotation and torsional motion leading to a breakdown 
of torsional alignment [12 . 

So far only very few experimental papers have investi- 
gated the possibility to control and do real-time mon- 
itoring of the torsion in axially chiral molecules. In 
Refs. [MlIT^ the stereogenic axis of a 3,5-difluoro-3',5 '- 
dibromo-biphenyl molecule was held fixed-in-space by 
1-dimensional (ID) adiabatic alignment with a linearly 
polarized 9 ns long laser pulse. In equilibrium the an- 
gle between the F-phenyl plane and the Br-phenyl plane 
is ±39°. Upon irradiation with a 700 femtosecond (fs) 
long, intense (but non- fragmenting), nonresonant pump 
pulse, torsion was induced accompanied by overall ro- 
tation around the fixed stereogenic axis. The torsional 
and rotational dynamics, measured by inducing Coulomb 
explosion with a short intense delayed probe pulse and 
recording of the emission direction of recoiling Br+ and 
F+ ions through ion imaging, showed that the ampli- 
tude and period of torsion was 0.6° and 1 picosecond 
(ps). While the period was in qualitative agreement with 
the theoretical model prediction of 1.2 ps, the amplitude 
predicted by theory was larger, 2.45°. The theory ratio- 
nalized that stimulated Raman transitions, driven by the 
pump pulse, creates a wave packet of torsional eigenstates 
in the electronic ground state leading to the torsion ob- 
served. It is one of the aims of the present work to resolve 
the discrepancy between theory and experiment. 

The principle of the current work is similar to the pre- 
vious work, [m [15] but three key experimental factors 
were changed to achieve much better resolved peaks in 
the experimental ion images which allows for tracing of 
the torsional motion with unprecedented precision and, 
therefore, a quantitative rather than a qualitative com- 
parison to the calculated results. The three decisive fac- 
tors are: 1) A chemically much purer sample resulting 
in higher contrast ion images. 2) The molecules were 3D 
aligned [161 HZ] prior to the kick and the probe pulses, 
i.e., the entire molecule was fixed-in-space not just the 
stereogenic axis (ID alignment). This implies that the 
torsional motion stands out much more distinctly in the 
ion images. 3) The molecule employed here, 3,5-difluoro- 
3',5'-dibromo-4'-cyanobiphenyl, is almost identical to 
the molecule studied in Refs. [l4l[T5] except that it has a 
nitrile group attached to the end of the Br-phenyl. This 
gives the current molecule a dipole moment of 4.4 Debye 
[18] compared to 0.2 Debye |T8] for the molecule used 
previously. The large permanent dipole moment of the 



molecules made it possible to use an electrostatic deflec- 
tor to spatially separate molecules according to their ro- 
tational quantum state [191 [20]. Performing the separa- 
tion prior to the interaction with the laser pulses, allowed 
us to focus the laser pulses onto only the molecules with 
low rotational energy, improving the degree of prealign- 
ment [19ii20j. 

The paper is organized as follows. In Sec. [n| the 
experimental setup is described. In Sec. |III| the results 



are discussed and Sec. IV provides a comparison with 
theory. Section [Vj concludes. 



II. EXPERIMENTAL SETUP 

The experimental setup has recently been described in 
detail [21 and only a brief description of the most impor- 
tant features is given here. A cold molecular beam (~1K) 
is formed by heating less than fifty milligram of solid 
3,5-difluoro-3 ',5 '-dibromo-4 '-cyanobiphenyl in 90 bar of 
He to 170° C and expanding the mixture into vacuum 
through a pulsed Even-Lavie valve. This chemical is not 
commercially available and was synthesized specially for 
the experiment (see Appendix). Here it should be noted 
that in the previously reported experiment the molecular 
sample was slightly contaminated which can be seen di- 
rectly from close inspection of the NMR spectrum in Fig. 
4 (d) of [15 . In particular the NMR spectrum shows a 
small impurity (5%) of the biproduct where both phenyl 
rings have flourine substituents. This contamination is 
problematic since the vapor pressure of the biproduct is 
higher than the target molecule and, therefore, the ratio 
between impurity and target molecules will be higher in 
the molecular beam than the initial 5% in the powder 
sample. Furthermore, the second most polarizable axis 
of the impurity will be located halfway between the two 
phenyl rings and, therefore, contribute to significant blur- 
ring of the F+ ion images and, thus, reduce the capability 
to experimentally resolve the torsional motion. 

The beam of cold molecules pass through an elec- 
trostatic deflector that deflects the molecules according 
to their effective dipole moment. The effective dipole 
moment is determined by its rotational quantum state 
[20 . Following deflection, the molec ules ent er a velocity 
map imaging (VMI) spectrometer (Fig. 1) where they 
are irradiated by three pulsed laser beams. The first 
pulse (YAG pulse: A = 1064 nm, tfwhm = 10 ns, 
Iyag = 6 X 10^^ W/cm^) serves to adiabatically align the 
molecules. It is elliptically polarized with the major axis 
parallel to the spectrometer axis, the Z-axis (see Fig. 1). 



The second pulse (kick pulse: 800 nm, 200 fs, Ikick = 
2 X lO^^W/cm^) initiates torsion and it is linearly po- 
larized along the Y-axis. The third pulse (probe pulse: 
800 nm, 30 fs, Iprobe = 3 x 10^^ W/cm^) is used to char- 
acterize the spatial orientation of the two phenyl planes 
by Coulomb exploding the molecules and recording the 
recoiling F+ and Br+ fragments. It is linearly polarized 
either along the Z-axis or the Y-axis. By recording series 




FIG. 1. Schematic of the velocity map imaging spectrome- 
ter used to detect F^ and Br^ ions. For illustrative purposes, 
the regions on the detector screen where F^ (Br^) ions are 
expected to hit for perfectly aligned molecules, at equilib- 
rium, are marked in yellow (black). The images behind the 
detector are experimental data recorded 2 ps after the kick- 
pulse. The polarization state of the YAG (alignment) pulse 
(A = 1064 nm, tfwhm = 10 ns, Iyag = 6 x 10^^ W/cm^) and 
the kick pulse (800 nm, 200 fs, Ikick = 2 x lO^^W/cm^) are 
also displayed. 



of images for different kick-probe delays, t, the evolution 
of the torsional motion can be monitored. To ensure that 
the molecules probed are both aligned and torsionally ex- 
cited, the spotsize of the three laser pulses are adjusted 
such that the foci of the YAG and kick pulses are large 
{ujQ = 38 pm and ujq = 35 pm) compared to the probe 
pulse {ujQ = 26 pm). 



III. EXPERIMENTAL RESULTS 

A. Deflection 

The effect of the deflector on the molecular beam is 
shown in Fig. 2| by the vertical intensity profiles. They 
are obtained by recording the magnitude of an ion sig- 
nal, produced by the probe pulse only, as a function of 
the vertical position of the laser focus similar to what was 
previously reported (see, for instance, Refs. pT) IT9| |2T] ) . 
Here, we use the strongest peak in the ion time-of-fiight- 
spectrum which is the C+ signal. When the deflector is 
off, the molecular beam extends over '^1.5 mm, mainly 
determined by the diameter of the skimmer before the 
deflector. When the deflector is turned on, the molecular 
beam profile broadens and shifts upwards (positive Y- 
values). The molecules in the lowest rotational quantum 
states have the largest effective dipole moments and are, 
therefore, deflected the most, as shown in recent work on 
iodobenzene and benzonitrile molecules [19l [20] . In the 
measurements described below, experiments were con- 
ducted on quantum state-selected molecules simply by 
positioning the laser foci close to the upper cut-off re- 
gion in the 10 kV profile, indicated by the red arrow in 
(Y = 1.3 mm). Our motivation for using the de- 
part of the beam is that despite the density of 
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FIG. 2. Vertical profiles of the molecular beam measured by 
recording the C^ signal as a function of the vertical position 
of the probe beam focus. The experimental data are shown by 
black squares (deflector off, kV) and red circles (deflector on, 
10 kV). The red arrow at Y = 1.3 mm indicates the position 
of the laser foci used to acquire ion images of the deflected 
molecules. 



B. Adiabatic prealignment 

The 3,5-difluoro-3 ' ,5 '-dibromo-4 '-cyanobiphenyl 

(CisHsFsBrsN, to be denoted by DFDBrCNBph in the 
following) is a chiral asymmetric top molecule, and the 
molecular composition in the supersonic expansion is 
a racemate, i.e., consisting of 50% of each of the two 
enantiomers. 

Three-dimensional alignment of the molecules is 
achieved using an elliptically polarized YAG pulse which 
confines both the most polarizable axis (MPA) and the 
second most polarizable axis (SMPA) of the molecule 
[li[l7l[22l|23], respectively. In DFDBrCNBph the MPA 
is located along the stereogenic axis (C-C b ond axis con- 
necting the two phenyl rings - see Fig. 1 ) whereas the 



SMPA is perpendicular to this at an angle of 10° from the 
phenyl-plane with the nitrile and bromine substituents 
towards the other phenyl ring (see Fig. 4). The effect 



of applying the YAG pulse, having its major and minor 
polarization axis along the Z-axis and the Y-axis, respec- 
ti vely, ca n be seen in panels (a) and (b) at t = — 0.33ps 



Fig. 3, showing F+ and Br+ ion images from Coulomb 



Fig. 2 
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molecules is lower, the alignment is higher pT f [191 120]. 



exploding the molecules at the peak of the YAG pulse. 
In both images two distinct features are observed. First, 
the center of the detector is almost void of ions (a minor 
region has been cut in the F+ images to exclude a noise 
signal from residual water molecules in the chamber). 
This observation is only compatible with the stereogenic 
axis (C-C bond axis) being confined perpendicular to the 
detector plane, i.e., along the Z-axis. Secondly, it is seen 
that both the Br+ and the F+ ions localize around the 
minor polarization axis of the YAG pulse (the Y-axis), 
the confinement being more pronounced for the Br+ ions. 
We interpret this as the SMPA of the molecule being 



aligned along the Y-axis. The more spread-out signal of 
the F+ ions results from the fact that the F-phenyl ring, 
at equilibrium, is offset 28° from the SMPA whereas the 
Br-phenyl ring is only offset by 10°. Thus, the joint ob- 
servations of the F+ and Br+ images demonstrate that 
the molecules are 3D aligned as expected. 



C. Time-dependent Torsion 

To investigate the torsional motion induced by the kick 
pulse we recorded series of both F+ and Br+ images at 
different kick-probe delays. The probe pulse was polar- 
ized perpendicular to the detector plane to ensure a de- 
tection efficiency that is independent of the orientation 
of the F- and Br-phenyl-rings. The results are shown in 
panels A and B of |Fig. 3[ The first striking observation is 
that from 0.67 ps to 4.0 ps the F+ images develop a pro- 
nounced 4-peak structure consisting of a pair in the top 
and a pair in the bottom part of the images. This shows 
that the kick pulse sharpens the alignment of the SMPA 
in agreement with a recent study on naphthalene where 
it was demonstrated that the combined action of 3D adi- 
abatic prealignment and a short kick pulse could signif- 
icantly improve the alignment of the molecular plane in 
a time interval shortly after the kick pulse [24 . The sec- 
ond observation is that the angular separation of the two 
regions in each F+ pair oscillates as a function of time. 
At t = 1.33 ps the angular separation reaches a local 
minimum, at t= 2.0 ps a local maximum, at t = 2.67 ps 
a local minimum etc. Turning to the Br+ ion images it 
is seen that the confinement along the Y-axis (vertical) 
improves. Unlike the F+ images a clear 4-peak structure 
is not observed. We ascribe this to the fact that the Br- 
phenyl ring is much closer to the SMPA, i.e., Br+ ions 
from molecules with the Br-phenyl ring localized on each 
side of the SMPA will overlap and prevent a distinct 4- 
peak structure. It is, however, observed that the width 
of the Br+ oscillates as a function of time in phase with 
that observed in the F+ images. These joint observations 
show that the dihedral angle between the Br-phenyl ring 
and the F-phenyl ring changes as a function of time. In 
other words, the torsional motion is directly imprinted 
on the ion distributions recorded. 

The contrast can be improved significantly by record- 
ing the F+ and Br+ images with the probe pulse polarized 
parallel to the kick pulse polarization. These data are 
displayed in Fig. 3, panels (c) and (d). The evolution of 
both ion species is the same as that observed in rows (a) 
and (b), but the contrast is significantly higher. In par- 
ticular, the 4-peak structure in the F+ images is very dis- 
tinct with complete separation between the two regions 
in each pair and the oscillations in the angular width of 
the Br+ images stand out very clearly. Our qualitative 
understanding of the improved contrast is that the probe 
pulse preferentially ionizes molecules with their SMPA 
close to its polarization axis. In other words, the probe 
pulse selects the molecules with the strongest alignment 



of the SMPA. 

The quantitative analysis is performed by fitting the 
angular distributions, extracted from the ion images, 
with a sum of four Gaussian functions. From these fits, 
the average angles, (0Br+)7 of the Br-phenyl ring, with re- 
spect to the kick pulse polarization, can be extracted di- 
rectly, from the peak positions of the Gaussian functions. 
Similarly, the average angle, (^f+)5 of the F-phenyl-ring 
can be found from the fit of the F+ angular distributions. 

Figure [4] displays the average angle between the F- 
phenyl-rings (black squares) and the kick pulse polar- 
ization as a function of kick-probe time delay for the 
time interval where a clear four-peak-structure can be 
identified in the angular distributions. The dashed black 
line splined to the experimental measurements indicates 
that the oscillations observed in the images of the F- 
phenyl ring have an amplitude of ~2.5°with respect to 
the kick pulse polarization direction and a period of 
1.25 ps. Closer examination of the oscillating trace re- 
veals a small overall increase towards larger (0f+) which 
shows that the plane of the F-phenyl-ring slowly moves 
away from the kick pulse polarization. Also the average 
angle between the Br-phenyl-rings (red circles) and the 
kick pulse polarization is shown. In the time interval 
1-4 ps the uncertainties on the measurements are small 
and allows for determination of the period (1.25 ps) and 
amplitude (~0.5°) of the motion. For the Br-phenyl-ring, 
we observe a small decrease in the (^Br+ ) with increasing 
time delay. The slight decrease in ((/)Br+ ) and correspond- 
ing increase in ((/>f+), are fully compatible with an overall 
rotation of the molecule about the stereogenic axis. 

The average values, ((/>f+) and ((/>Br+)7 can be used 
to calculate the instantaneous dihedral angle, ^d, of the 
molecule by simple addition (see Fig. 4 (b)). By retriev- 
ing these average values at every kick-probe time delay it 
is possible to directly monitor the evolution of the dihe- 
dral angle of the molecule. The results for the dihedral 
angle ((/)d, blue triangles) are displayed in Fig. 4 (a), with 



a dashed line included to help guide the eye and empha- 
size the torsional motion. The oscillation is seen to have 
an amplitude of 3° and a period of 1.25 ps. 

We also recorded F+ and Br+ ion images for times 
larger than 4 ps. Th e results, including the 3.0 ps data, 
are shown in Fig. 5 For t > 4.0 ps the F+ ion images 
loose their distinct four region structure and, instead, 
evolve into a circular shape at t > 8 ps. This prevents 
tracking the evolution of (^f+)- The images, however, 
confirm that the stereogenic axis of the molecule is still 
tightly confined along the major axis of the YAG po- 
larization ellipse, since the central region is absent from 
ions. The corresponding Br+ ion images show a simi- 
lar delocalization behavior, but the broadening is slower. 
This observation is consistent with theory as discussed in 
ISec. IV Bl 

To quantify the time-dependence of the delocalization 
in the F+ and Br+ ion images we determined for each 
image (cos^ 0^21^)7 where a2D is the angle between the 
projection of an ion velocity vector on the detector plane 
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FIG. 3. F^ (panel (a) and (c)) and Br^ (panel (b) and (d)) ion images obtained as a function of the kick-probe time delay 
indicated in ps above each column of panels. Panels (a) and (b) ((c) and (d)) were obtained using a linearly polarized probe 
pulse being perpendicular to (parallel with) the detector. The delay of the probe pulse (A = 800 nm, tfwhm = 30 fs, Iprobe 



= 3 X lO^'^W/cm^) with respect to the kick pulse (A 
numbers on the top of the vertical panels. 



800 nm, tfwhm = 200 fs, Ikick = 2 x 10 '^W/cm'') is given by the 



(^) SMPA 
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FIG. 4. (a) Time-dependence of the dihedral angle (blue tri- 
angles) following the kick-pulse, obtained by adding the orien- 
tation of the F-phenyl-ring (black squares) and the orientation 
of the Br-phenyl-ring with respect to the SMPA - see text for 
details, (b) Sketch of the molecular structure as seen in end 
view. The SMPA of the molecule and the extracted angles 
are marked on the drawing. 



F-phenyl-ring 
Br-phenyl-ring 




3 6 9 12 

Kick-probe delay (ps) 



15 



FIG. 6. Time dependence of the planar alignment of the in- 
dividual phenyl-rings. The (cos^a2D) values are calculated 
from rotated angular distributions to remove the effect of the 
torsion (see text for details). Black squares (red circles) cor- 
respond to the confinement of the F-phenyl-ring (Br-phenyl- 
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FIG. 5. F^ (a) and Br^ (b) ion images obtained at long 
{t > 3 ps) kick-probe time delays obtained using a probe pulse 
linearly polarized perpendicular to the detector. The delay 
of the probe pulse (A = 800 nm, tfwhm = 30 fs, Iprobe = 
3 X 10^^ W/cm^) with respect to the kick pulse (A = 800 nm, 
TFWHM = 200 fs, Ikick = 2 X lO^^W/cm^) is given by the 
numbers on the top of the vertical panels. 



and the Y-axis (See [Fig. 7 



This is implemented by ro- 
tating the peaks fitted to the angular distributions of F+ 
and Br+ ions, towards the Y-axis by the following proce- 
dure: If the peak is located in the first or third quadrant 
the F+ (Br+) image is rotated by — (^p+) (— (0Br+)) be- 
fore calculation of (cos^a2D), whereas if the peak is in 
the second or fourth quadrant the image is rotated by 
+ ((/>F+) (+(^Br+))- This procedure allows for calcula- 
tion of the best estimate of the spread from the mean by 
including all data. The (cos^Qf2D) values are displayed 
A clear difference in the delocalization dynam- 



m 



Fig. 6 



ics of the two phenyl planes is observed, with the light 



fluorine substituted phenyl ring approaching (cos^Qf2D) 
= 0.5 more rapidly than the heavier bromine substituted 
phenyl ring. This is fully consistent with semi-classical 



theoretical predictions presented in Sec. IV B 



D. Covariance and Autovariance Mapping 



Before discussing the theoretical results we introduce 
covariance mapping analysis [25 to the experimental 
data. This analysis strengthen our interpretation of tor- 
sional motion, based on the fluorine and bromine ion 



images in Fig. 3 Covariance mapping is a technique to 



reveal correlations, which in the high count rate regime 
would otherwise become blurred or completely lost due 
to the large contributions from uncorrelated events. The 
basic principle is to calculate the covariance, i.e., to ob- 
tain a crosscorrelation of the variance in the data [26 , to 
extract the deviation from the mean value for all the indi- 
vidual frames of the ion image, pertaining to each (probe) 
laser shot. This is fundamentally different from coinci- 
dence measurements in, e.g., COLTRIMS and reaction- 
microscopes |27l |28] where the count rate is restricted to 
less than one event per laser shot, since the covariance 
map is produced from the correlations of the variance, 
which per definition has no restrictions on the number of 
events obtained per laser shot. 

The covariance {C{x,y)) between two observables 
X{x) and Y{y) is defined as the mean of the product 
between the observables {X{x)Y{y)) subtracted by the 
product of the means {X (x)) {Y (y)) : 



C{x,y) ={{X{x) - (XixMYiy) - (F(y)))) 
={X{x)Y{y)) - {X{x)){Y{y)) 



(1) 



Here X and Y correspond to the observed signals, with 
X and y being the variables investigated for correlation. 

Although many different correlation techniques have 
been explored previously and found many different ap- 
plications within natural science [28H3T] . to our knowl- 
edge no reports on using covariance mapping to extract 
angular correlation information from 2D ion images has 
been reported. We suggest to do this by replacing X 
and Y with O and S, corresponding to the angular dis- 
tributions of the ion signal, and x and y with 6 and ^ 
being th e ion detection angle with respect to the Y-axis 
(Fig. 1). Applying this version of covariance mapping 



to the impact coordinates of the ions in the 2D velocity 
map images, allows for extraction of the angular correla- 
tion between ejected fragments. In our experiment only 
one ion species, e.g., F+, is detected per laser shot and 
thus O = S. As a consequence Eq. M may be written 




FIG. 7. The autovariance mapping principle illustrated for 
F^ fragments produced by coulomb explosion of laser aligned 
DFDBrCNBph. From the autovariance map seen on the right 
it is possible to extract correlation information between the 
ejection angle of the F^ ions. The blue arrows indicate the 
steps in the data acquisition and data processing leading to 
the production of autovariance maps. 



as 



c{0u02) =(6(^1)6(^2)) - (e(^i))(6(^2)) 
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Since this corresponds to investigating correlations of a 
signal with itself this approach is equivalent to obtaining 
an autocorrelation of the variance in the data and so the 
2D maps depicting the correlation should be termed au- 
tovariance maps. The principle of autovariance mapping 
is illustrated in Fig. 7\ Following each laser shot the ions 
focused onto the detector are recorded by a CCD camera 
after which an image analysis program saves the impact 
coordinates, illustrated by the yellow bars in the image 
analysis box. The experiment is repeated N times, typ- 
ically N = 10.000, giving the averaged ion images (such 
as those shown in Fig. 3) with the averaged angular dis- 
tributions (Fig. 7 dashed black line). The autovariance 
C {61^62) can now be determined by applying Eq. ([2| 
to all N frames of the ion image, and displayed as a 2D 
map of the ejection angles. In this map the two axes cor- 
respond to the angular distribution from the averaged 
image, with the positive Y-axis denoting 0°. 

Autovariance maps obtained through this procedure 
will tend to produce a strong autocorrelation line along 
the diagonal of the map, since the observation of an ion in 
frame 6^ at the angle 61 is naturally also observed in the 
same frame, for the identical observable 6^, at O2 = Oi. 
This autocorrelation trace will also act as an axis of sym- 
metry for the correlation, since C(^i,^2) = C{92-,0i)^ 
which is in fact equivalent to the square of the r.m.s. 
standard deviation by definition. The positive areas of 
the autovariance maps, which fall outside the autocor- 
relation line (6^1 7^ ^2), indicate the correlation between 



ions at 6i and O2 , showing an increased likelihood of ob- 
serving an ion at O2 provided that an ion is detected at 

Figure [S] shows the angular autovariance maps ob- 
tained from both F+ and Br+ ion images at two different 
kick-probe time delays. In these maps the dominant au- 
tocorrelation signal along the diagonal has been set to 
zero to increase the contrast of the interesting but less 
intense correlations in the autovariance map. Also these 
data sets have been binned in steps of 5° to smoothen out 
small fluctuations. A striking feature of all autovariance 
maps is that the observed correlations are restricted to 
areas where the angle between the ejected ions is close 
to 180°. This is seen, e.g., in the F+ map at 2 ps, where 
the signal at ~150°(~210°) is correlated to the region at 
~-30°('^30°). In other words, when an F+ ion is ejected 
upwards at an angle of 150° there is a high probability 
that another F+ ion is ejected at -30°, i.e., just in the 
opposite direction. This corroborates our interpretation 
of the two ion species as direct observables of the orienta- 
tion of each of the phenyl rings. In addition, the observa- 
tion of four prominent positive areas in the autovariance 
map at t = 2 ps confirms the interpretation of the 4-peak 
structure in the F+ ion images (Fig. 3). Similar consid- 
erations hold for the Br+ autovariance map at t = 2 ps. 
Notably the appearance of a 4-peak structure with posi- 
tive covariance shows even clearer than the averaged ion 



arb. units 



image on Fig. 3, that the Br-phenyl plane has a small 
angular offset with respect to the kick pulse polarization. 
At t = 10 ps there is no longer a 4-peak structure but the 
distinct oblique lines show unambiguously that emission 
of both F+ and Br+ ions still occurs in a pairs with an 
upward and a downward ion. The extension of the co- 
variance principle to obtain angular correlation between 
ions allows for the extraction of additional information 
from the 2D images obtained using VMI spectrometers. 
As shown here, a single detector setup, can be used to 
obtain autovariance maps between ions of identical mass- 
to-charge ratios, which can be used to both substantiate 
and clarify the interpretation of the recorded ion images. 
This can be extended to obtain the correlation between 
ions with different mass-to-charge ratios by pulsing the 
front of the MCP twice per laser shot and running the 
CCD camera in burst mode. In addition, we note that 
the technique can also be extended to obtain radial auto- 
variance maps, for fragmentation processes with multiple 
energy channels, or combined to encompass both angular 
and radial correlation of the 2D ion images. 



IV. COMPARISON WITH THEORY 



The main findings of |Sec. III| were that the dihedral 
angle oscillates with a period of approximately 1.25 ps, 
and an amplitude of 3°. These findings were based on 
an analysis of the data assuming that the two ion species 
move out to the detector in the plane of the (substituted) 
phenyl ring from which they originated. The covariance 
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FIG. 8. Autovariance maps for F^ (upper panels) and Br^ 
ions (lower panels) at kick-probe time delays of 2 ps (left 
column) and 10 ps (right column). To smoothen out small 
fluctuations, the data has been binned in steps of 5 degrees. 
Also the large signal, due to autocorrelation, along the diago- 
nal of the maps, has been set to zero to reveal the less intense 
correlations in the autovariance map. 



and autovariance analysis introduced validated this as- 
sumption and we may build the theory on this fact. 



A. Early time dynamics: Torsional motion 

The purpose of this section is to provide a theoreti- 
cal foundation for understanding the experimentally ob- 

^ Atomic units (\e\ = 

h = n 



served ion images shown in Fig. 3 

e = ao = 1) are used throughout. A useful starting 
point is to discuss the widths of the ion images. These 
may conveniently be quantified in terms of the variance 



f)-(0i)2, i = Br+,F^ 



(3) 



For the sake of interpretation, we rewrite these two 
angles in terms of the dihedral angle, (j)^^ and the overall 



Fig. 4 (b)]: 




</'Br+ =^ + V<l^d, 


(4) 


</.p+ = $ - (1 - r?)0d, 


(5) 


._ -^F^ 


ff{\ 



Consequently, we may express the variance of the ob- 



served ions as 



2 I 2 2 



2r,(($(/.d) - (*)(<^d)), (7) 

4 + (1 - vfal^ 

2(l-,?)(($(/.d)-($)(<^d)). (8) 



4,^ = 4 + (1 - ^)'<^; 



Previous investigations show that the variance of the 
dihedral angle is small (See Figs. 3 and 9 in Ref. [15 ). 
It is therefore obvious to ask if cr| :^ a^^. As follows 
from Eqs. (JTl) and ([s]) this can only be the case if the 
widths of the ion images are comparable. Inspection of 
Figs, [s] and [g] confirms that for long times {t > 3 ps) this 
is not true. For shorter times, however, the widths of 
the two ion images are close (see Fig. 3). Further, we 
have previously argued for separability of the ^ and (j)^ 
coordinates in this limit [15 meaning that the last term 
in Eqs. ^ and ^ vanishes completely. In this regime 
$ is an adiabatic parameter and we simply discuss the 
short-term dynamics based on calculations with <I> fixed 
meaning that only torsion takes place. 

The torsional motion is then dictated by the time- 
dependent Schrodinger equation 

idt^i^; ^d, t) = [Td + VtoriM + ^kick(^; (/>d, t)] ^($; (/)d, t), 

(9) 
with 



Td = 






(10) 



being the rotational kinetic energy due to torsion, Ftor 
the torsional potential and Vkick the polarizability inter- 
action energy between the kick pulse and the molecule: 

^kick(^,0d,^) = -^i^o(^)[Q^xx(^d)cos^(^ + ^^d) 

+ %y(</>d)sin^(*^ + ^</>d) 
- 2axy((^d) cos(<I> + 7?(^d) sin(<l> + rjcj)^)]. 

(11) 

We solve this equation for DFDBrBPh using a close cou- 
pling method with a value of ^ corresponding to the 
SMPA and kick pulse polarization aligned. As initial 
state ^(to), we use a state localized in the torsional well 
at 39°. The details have been given in Ref. [15]. 

Figure |9] shows the result of the calculation. The tor- 
sional dynamics consist of small, periodic oscillations. As 
argued previously [TF, these oscillations stem from the 
fact that the kick pulse leaves DFDBrBPh in a coherent 
superposition of the ground state and first excited state of 
the torsional potential — essentially simple harmonic os- 
cillator states separated by 3.42 meV. The molecule used 
in the experiment is DFDBrCNBPh rather than DFD- 
BrBPh. The torsional potential differs, however, only 
slightly from that used previously for DFDBrBPh. Since 
the kick pulse only populates the first few torsional states 
where the potential is very close to harmonic, the main 
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FIG. 9. Comparison: Blue triangles - experimental results, 
dashed blue line - theoretical calculation on DFDBrBPh 
shifted by -0.9° to match the dihedral angle of DFDBrC- 
NBPh at equilibrium. 



effect to take into account is the shift of the minima of 
the torsional potential from one molecule to the other 
(—0.9°). It follows that we may simply shift the pre- 
viously calculated theoretical curve by —0.9° to obtain 
a comparison with the experimental results. As seen in 



Fig. 9 this leads to an impressive agreement with the 
experimental results. 



B. Long time behavior: Dephasing of the phenyl 
rings 

For longer time scales {t > 3 ps), the semi-classical 
model is no longer attractive. The reasons are that the 
widths of the two ion images (Figs, p^ and p| start to 
deviate from one another meaning that the variance is 
no longer dominated by $ alone and, hence, that the 
separability in ^ and 0d breaks down [12] |T5] . 

To provide a qualitative analysis in this long-time 
regime, we exploit a different theoretical limit, where the 
phenyl-rings are treated as rigid rotors. In this case, each 
of the rings i = Br^,F^ can be described by a wave 
packet 



*i( 






27r 



eiJ.^ig- 



(12) 



where Ji labels a rotational state and where cj. is a fixed 
constant after the end of kick pulse. 

To estimate on what time scale, ri, such a wave packet 
dephases, we look at the time evolution of the coherence 
between two neighboring levels Ji and Ji + 1. It immedi- 
ately follows from (12) that this coherence beats like 



2Ji + l 
2/i 



(13) 



Consequently, the two rings dephase at a ratio 
rp+ /p+ 2 Jb^+ + 1 



^Br+ 



/Br+ 



2Ji 



F+ 



(14) 



The value of Ji depends on the rotational kinetic energy, 
£^kin,i, of each ring. Due to the a times higher polarizabil- 
ity of the Br-phenyl-ring compared to the F-phenyl-ring, 
we assume that £^kin,Br = ax £^kin,F with £^kin,i = J\ /"^h- 
This result is achieved by integrating the torque of a clas- 
sic rigid rotor from 0° to 90° to obtain an energy esti- 
mate. We then immediately have 



^Br+ 




(15) 



Using this in Eq. ( 14 ) yields 




(16) 



With a value of 0.21 for the ratio of the moments of 
inertia and a = 2 [18] we obtain r = 0.65 in qualitative 
agreement with the dephasing ratio of 0.68 obtained from 
the data presented in |Fig. 6| by comparing the full width 
at half maximum values of the Gaussian fits to the ion 
pictures. 



V. CONCLUSIONS 

In conclusion we have shown that it is possible to in- 
duce torsion in a substituted biphenyl molecule with a 
nonresonant kick pulse and image the motion in real-time 
with an intense delayed probe pulse. The measurements 
show a distinct torsional motion with an amplitude of 3° 
and a period of 1.25 ps for the first 4 ps after the kick 
pulse, in excellent agreement with results from our theo- 
retical model. At longer times delocalization of the two 
phenyl rings of the molecule blurs torsion. 

An important prerequisite for the experiment is the 
ability to keep the molecules fixed-in-space during the 
time-resolved experiment, practically obtained by 3D 
aligning them in the adiabatic limit. This places both 
benzene rings of the molecule perpendicular to a 2D 
imaging detector. Upon Coulomb explosion by the probe 
laser Br+ ions from one benzene ring imprints an ion im- 
age on the detector screen that uniquely identifies the 
orientation of this ring. Likewise, F+ ion images identify 
the orientation of the second ring. The difference between 
the orientation of the two rings determines the dihedral 
angle, i.e., the normal vibrational coordinate character- 
izing torsion. Such a procedure would not be possible for 
samples of randomly oriented molecules. 

The paper also showed that it is possible to apply co- 
variance mapping analysis of ion images. This revealed 
clear correlations between the emission direction of ions 
strengthening our interpretation of Coulomb explosion 



as a direct and unique observable of molecular orienta- 
tion. Covariance and autovariance methods applied to 
ion imaging have a potential to provide additional infor- 
mation beyond what is normally extracted from averaged 
quantities. 
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APPENDIX 

Synthesis of 
3,5-difluoro-3 ',5 '-dibromo-4 '-cyanobiphenyl 

3,5-Difluoro-4'-aminobiphenyl:[32l To a sample vial 
was added 4-bromoaniline (172.0 mg, 1.0 mmol), 3,5- 
difluorophenyl boronic acid (189.5 mg, 1.2 mmol), tetra- 
butylammonium bromide (32.4 mg, 0.1 mmol), potas- 
sium carbonate (276.4 mg, 2 mmol) and water (2 mL). To 
this was added 100 pL of a PdClsiEDTAiNasCOa stock 
solution prepared by mixing PdCl2 (17.7 mg, 0.1 mmol), 
EDTA (37.2 mg, 0.1 mmol) and Na2C03 (21.2 mg, 0.2 
mmol) in water (1 mL). The vial which was fittet with a 
Teflon sealed screwcap and heated to 100° C for 45 min 
after which the reaction was cooled to room tempera- 
ture. The mixture was diluted with Et20 (25 mL) and 
CH2CI2 (10 mL). The organic phase was washed with wa- 
ter (2x15 mL), brine (15 mL) and dried using MgS04. 
The solvents were removed under reduced pressure and 
the crude reaction mixture was purified by column chro- 
matography eluting with CH2CI2 to afford the title com- 
pound as a colorless solid (170.5 mg, 83%). ^H NMR 
(400 MHz, CDCI3) 6 (ppm) 7.37 (d, 2H, J = 8.2 Hz), 
7.07-7.00 (m, 2H), 6.74 (d, 2H, J = 8.2 Hz), 6.69 (tt, IH, 
J = 8.8, 2.3 Hz), 3.79 (bs, 2H). ^^C NMR (400 MHz, 
CDCI3) 5 (ppm) 163.6 (dd, J = 245.0, 13.2 Hz), 147.3, 
144.8 (t, J = 9.6 Hz), 129.1 (t, J = 2.5 Hz), 128.2, 115.6, 
109.1 (dd, J = 18.4, 6.9 Hz), 101.5 (t, J = 25.4 Hz). ^^F 
NMR (376 MHz , CDCI3) 5 (ppm) -110.8 (t, J = 8.8 
Hz). HRMS. C12H9F2N [M+H+] Calculated: 206.0781. 
Found: 206.0780. 

3,5-Difluoro-3 ' ,5 '-dibromo-4 '-aminobiphenyl: ^ 3,5- 
Difluoro-4 '-aminobiphenyl (165.0 mg, 0.8 mmol) was dis- 
solved in glacial acetic acid (2.5 mL). Bromine (82iiL, 
255.7 mg, 1.6 mmol) dissolved in glacial acetic acid (3.8 



10 



mL) was slowly added. The reaction was left stirring 
for 1 hour. Then water (5 mL) and CH2CI2 (30 mL) 
was added. Sat. Na2C03 was added until the reaction 
mixture became alkaline. The organic phase was washed 
with sodium thiosulfate (10 mL), water (10 mL), brine 
(10 mL) and dried using MgS04. The organic solvent 
was removed under reduced pressure affording the title 
compound as a pale brown solid (271 mg, 93%). This 
was used without further purification. ^H NMR (400 
MHz, CDC13) 6 (ppm). 7.59 (s, 2H), 7.01-6.96 (m, 2H), 
6.74 (tt, IH, J = 8.84, 2.3 Hz) 4.69 (bs, 2H). ^^C NMR 
(400 MHz, CDCI3) 6 (ppm) 163.6 (dd, J = 246.6, 13.1 
Hz), 130.4 (t, J = 2.6 Hz), 130.3, 109.29 (dd, J = 18.6, 
7.2 Hz), 109.20, 102.6 (t, J = 25.2 Hz). ^^F NMR (376 
MHz , CDCI3) S (ppm) -109.8 (t, J =8.5 Hz). GCMS. 
Ci2H7F2NBr2 Calculated: 362.99 Found: 363 (100), 203 
(36), 175 (24), 101 (14). 

3,5-difiuoro-3 ' ,5 '-dibromo-4 '-cyanobiphenyl: Re- 

action mixture 1: 3,5-Difiuoro-3',5'-dibromo-4'- 
aminobiphenyl (270.0 mg, 0.74 mmol) was dissolved 
in water (1.0 mL) and glacial acetic acid (3.7 mL). 
Concentrated sulfuric acid (1.0 mL) was added and the 
reaction mixture was cooled to 10° C. Sodium nitrite 
(56.5 mg, 0.82 mmol) dissolved in a minimum of water 



was slowly added and the reaction was left stirring for 
30 min at which point the reaction mixture had become 
clear. Reaction mixture 2: In a second round bottom 
fiask CuS04*5H2 (223 mg, 0.89 mmol) is carefully 
mixed in water (1 mL) with KCN (242 mg, 3.72 mmol) 
dissolved in water (ImL) keeping the temperature below 
20° C cooling with ice. To reaction mixture 2 is slowly 
added reaction mixture 1 while keeping the reaction 
basic by adding Na2C03 (sat.). After the addition is 
complete the reaction is left stirring 1 hour at room 
temperature. CH2CI2 (30 mL) and Et20 (50 mL) are 
added and the organic phase was washed with water 
(2x20 mL), brine (20 mL) and dried using MgS04. 
The solvents were removed under reduced pressure and 
the crude reaction mixture was purified by column 
chromatography eluting with pentane:CH2Cl2 3:1 to 
afford the title compound as a colorless solid (141.8 mg, 
50%). ^H NMR (400 MHz, CDCI3) S (ppm) 7.8 (s, 2H), 
7.11-7.05 (m, 2H), 6.92 (tt, IH, J = 8.6, 2.3 Hz). ^^C 
NMR (400 MHz, CDCI3) S (ppm) 163.7 (249.1, 12.7 
Hz), 145.2, 140.0 (t, J = 9.4 Hz), 130.5, 127.5, 118.4, 
115.9, 110.7 (dd, J = 18.9, 7.7 Hz), 105.3 (t, J = 25.0 
Hz). i^F NMR (376 MHz , CDCI3) S (ppm) -107.7 (t, 
J = 8.0 Hz). GCMS. Ci3H5F2NBr2 Calculated:372.99. 
Found: 373 (100), 213 (68), 186 (15), 106 (15). 
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